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Abstract 

Roughly speaking, a near-best (abbr. NB) quasi-interpolant (abbr. QI) is an approx- 
imation operator of the form Q a f = Yla£A-^-a{f)B a where the -B a 's are B-splines and 
the A a (/)'s are linear discrete or integral forms acting on the given function /. These 
forms depend on a finite number of coefficients which are the components of vectors a a 
for a G A. The index a refers to this sequence of vectors. In order that Q a p = p for 
all polynomials p belonging to some subspace included in the space of splines generated 
by the B^s, each vector a a must lie in an affine subspace V a , i.e. satisfy some linear 
constraints. However there remain some degrees of freedom which are used to minimize 
||a a ||i for each a & A. It is easy to prove that max{||a a ||i; a G A} is an upper bound of 
|| Q a || oo : thus, instead of minimizing the infinite norm of Q a , which is a difficult problem, 
we minimize an upper bound of this norm, which is much easier to do. Moreover, the 
latter problem has always at least one solution, which is associated with a NB QI. In the 
first part of the paper, we give a survey on NB univariate or bivariate spline QIs defined 
on uniform or non-uniform partitions and already studied by the author and coworkers. 
In the second part, we give some new results, mainly on univariate and bivariate integral 
QIs on non-uniform partitions: in that case, NB QIs are more difficult to characterize 
and the optimal properties strongly depend on the geometry of the partition. Therefore 
we have restricted our study to QIs having interesting shape properties and/or infinite 
norms uniformly bounded independently of the partition. 



2 



P. Sablonniere 



1. Introduction and notations 

The spline quasi-interpolants (abbr. QIs) considered in this paper have the following 
general form: 

A denotes a finite or infinite set of indices. B a is a B-spline with support S a , defined 
on a uniform or nonuniform partition. The index a of Q a refers to a family of vectors 
{a a ,a G A} which is described below. The coefficients A Q (/) are discrete or integral 
functionals of the following types: 




where T a is a finite set of indices nearby a. In the discrete coefficient functionals A a (f), 
the points x 1 lie in S 7 and the associated QIs are called discrete quasi-interpolants 
(abbr. dQIs). Those with integral coefficient functionals A Q (/) are called integral quasi- 
interpolants (abbr. iQIs): the weight function in the integral is a B-spline B* (which 
can be different from 5 7 ) satisfying J E B* = 1 for all indices 7. For a£i, the vectors 
a a = (a a (7), 7 G T a ) are determined in the following way: 

(i) let P be a space of polynomials included in the space of splines generated by the 
family {B a , a G A}, then we impose that Q a be exact on P, i.e. Q a p = p for all p G P. 
This is equivalent to say that, for each a G A, the vector a a belongs to some affine 
subspace V a of where n a = card(T a ) is large enough. 

(ii) in general there remain some undetermined components of a a . Then we minimize 
the ^-norm of this vector: using standard results in optimization, it is easy to prove that, 
for each a E A, there exists at least one vector a a G V a solution to this problem. 

The infinite norm of Q a is bounded above by v{Qa) = rnax a€j 4 v{a a ) where v(a a ) = 
||a a ||i- Thus, instead of minimizing this norm, we minimize z/(a Q ) for all a G A, the 
vector a a satisfying the constraints a a G V a , which is much easier to do. When a* is 
a solution of this minimization problem, we say that the corresponding QI Q* = Q a * is 
a near-best quasi-interpolant (abbr. NBQI). In the first part of the paper, we present 
various examples of NBQIs which have been studied in PP|22] and in the reports |2]-|5- 
In the second part, we present new results on some families of univariate and bivariate 
iQIs: the proofs are somewhat sketched and will be detailed elsewhere |S]|2BJi26j. 
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2. Univariate Near-best QIs on uniform partitions 

Butzer et al. [10] [33] have completed the results of Schoenberg |12] about the general 
expressions of spline dQIs of maximal approximation order (i.e. exact on polynomials 
having the same degree as the underlying spline) on the real line with Z as sequence of 
knots. Another technique, based on central difference operators, has been given by the 
author in p£fj . 

2.1. Near-best spline dQIs. Let M2 m (x) denote the centered cardinal B-spline. Con- 
sider the family of spline dQIs of order 2m depending on n + 1 arbitrary parameters 
a = (do, cl\i ■ ■ ■ , a, n ), n > m: 

QJ = ^A/(z)M 2m (x-z), 

with coefficient functionals 

n 

A/(i) = a f(i) + % (/(* + J) + /(< " •?')) ■ 

3=1 

Setting u(a) = \a \ + YTj=x \ a j\i then we have ||Q a ||oo < K )- By imposing that Q a be 
exact on the space IT of polynomials of degree at most r, with < r < 2m — 1, we obtain 
a set of linear constraints: a E V r G K n+1 . We say that Q* = Q a * is a near best dQI if 

v(a*) = min{z/(a); a G K-}. 

There is existence, but in general not unicity, of solutions (see [2~2~]). 

Example: cubic splines [22]- There is a unique optimal solution for r = 3 and n > 2: 

= 1 + 3^' a ™ = ~6^' ^* = ° f ° T 1 " J " n ~ L 
Moreover, for all n > 4, ||<5*||oo < 1 + 3^2- Here are the first values of ||<5*||oo & K a *)> 
n = 1 : 1.222 & 1.666; n = 2 : 1.139 & 1.166; n = 3 : 1.074 & 1.074. 

2.2. Near-best spline iQIs. A similar study can be done for integral spline QIs. We re- 
fer to 2\ [22] and we only give an example given in these papers. Setting a = (a , ai, . . . , a n ), 
n > m, Afi(x) = M 2m (x-i) and (/, M») = / /Mi, we consider Q a / = J2i& A f( l ) M i with 
coefficient functionals 

n 

A/(z) = a (/, Mi) + ^ aj ((/, Mi_,-> + (/, M i+j )) . 
3=1 

As in section 4.2, we have ||Q a ||oo < v ifl) an d we say that Q* = Q a * is a near best iQI if 
v{a*) = min{z/(a); a G V r }. There is existence, but in general not unicity, of solutions. 
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Example: cubic splines [22] • There is a unique optimal solution for r = 3 and n > 2: 

2 1 
°° = 1 + 3^ ' °™ = ~3^' Gj * = ° /0r 1 " j ~ U ~ L 
Moreover, for all n > 4, ||<3*||oo < 1 + 3^2- Here are the first values of ||Q*||oo & v(a*); 
n = 1 : 1.5278 & 2.333; n = 2 : 1.2778 & 1.333; n = 3 : 1.1481 & 1.1482. 



3. Bivariate Near-best QIs on uniform partitions 

3.1. A general construction of dQIs. Let <p be some bivariate B-spline on one of the 
two classical three or four direction meshes of the plane (e.g. box-splines or H-splines, 
see 00 00111 EH El citeSab6|H]). Let £ = supp(ip) and £* = £ n Z 2 . Let a be 
the hexagonal (or lozenge=rhombus) sequence formed by the values {(p(i),i G X*}. The 
associated central difference operator D is an isomorphism of P(<£>), the maximal subspace 
of "complete " polynomials in the space of splines S(ip) generated by the integer translates 
of the B-spline ip (see e.g. 00). Computing the expansion of a in some basis of the 
space of hexagonal (or lozenge) sequences amounts to expand D in some basis of central 
difference operators. Then, computing the formal inverse T>~ x allows to define the dQI 

Qf = V- l f{k) V {- - k) 

k£Z 2 

which is exact on P(<£>) . Let us now give two examples of NB dQIs which are detailed in 
[22J. The definition of these operators is quite the same as in section 2.1. 

3.2. Near-best spline dQIs on a three direction mesh. For example, let (p be the 
C 2 quartic box-spline with support S = H2 where H s denotes the regular hexagon with 
edges of length s > 1, centered at the origin and H* = if s flZ 2 . The near-best dQIs, which 
are exact on have coefficient functionals with supports consisting of the center and of 
the 6 vertices of H*, s > 1. The coefficients of values of / at those points are respectively 
l + o"^ an d — 12~2> therefore the infinite norm of the optimal dQIs Q* is bounded above 
by v* = 1 + 4r. Here are the first values of ||<2*||oo k v* s ; n = I : 1.34028 & 2; 
n = 2 : 1.22917 & 1.25; n = 3 : 1.10185 & 1.111. 

3.3. Near-best spline dQIs on a four direction mesh. For example, let <p be the 
C 1 quadratic box-spline. Let A s be the lozenge (rhombus) with edges of length s > 1, 
centered at the origin, and let A* = A s fl Z 2 . The near-best dQIs which are exact on 
n 2 have coefficient functionals with supports consisting of the center and the 4 vertices 
of A*, s > 1. The coefficients of values of / at those points are respectively 1 + ^ and 
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— gj2, therefore the infinite norm of the optimal dQIs Q* is bounded above by u* = 1 + 
Here are the first values of ||Q*||oo & n = 1 : 1.5 & 2; n = 2 : 1.25 & 1.25; 
n = 3 : 1.111 & 1.111. 

Examples with C 2 quartic box-splines are also given in [2*2"] . 



4. Univariate Near-best QIs on non-uniform partitions 

In this section, we consider dQIs or iQIs of degree m > 2 defined on a bounded interval 
/ = [a, b). Let E m = {— m + 1, . . . , 0} and J = {0, 1, . . . , m + n — 1}. Let T = {tjj z e 
_E m U J} be an arbitrary non-uniform increasing sequence of knots with multiple knots 
at t = a and t n = b, as usual. Let Bj be the B-spline with support [t,_ m , for 
j G J, and let e p (x) = x p for all p > 0. Setting 0j = — X^es V-M; ^ * s wen known that 

e i = EjejQjBj and e 2 = EjeJ 6 f B ^ with = T^hy £ (r ,.)si& > r<«VH-*i+«- We recall 
the following expansion [2*7] : 

Aj = 2 - = ^7 -T- (tj +r - tj+sf > 0. 

m 2 (m - 1) ^— ' 

(r,s)££2 i ,r<s 

More generally e r = Y^j£j@j B j wnere the flj s are proportional to symmetric functions 
of knots. The simplest dQI is the Schoenberg operator 

which is exact on Hi and is also shape preserving. Moreover, it satisfies 

S\e 2 - e 2 = E XjBi > 0. 

4.1. Uniformly bounded dQIs. Let us only give an example: we start from a family 
of differential QIs of degree m which are exact on LT2 (see also [3]): 

s 2f = E *? } (/)% A f (/) = M) - ^ 2 M)- 

On the other hand, ~D 2 f(6j) can be replaced on the space n 2 by the second order 
divided difference [6j-i,6j,0j + i]f, therefore the dQI defined by 

s *f = E^ftf)^ /4 2) (/) = m) - ..(>,,.(>.,■■: i 

is also exact on n 2 . Moreover, one can write 

/4 2) (/) = a ifi-i + b ifi + °ifi+i 
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with a % = -Ai/Afli-^Afli-j + A0i), Q = -Ai/A^A^ + A0«), and 
&i = 1 + Aj/ A6i-iA6i. So, according to the introduction 

H^Hoo < max(|aj| + + |cj|) < 1 + 2 max — . 

The following theorem [I] extends a result given for quadratic splines in |2~2lj [33j . 



Theorem 1. For any degree m, the dQIs S% are uniformly bounded independently of the 
partition. More specifically, if [r] denotes the floor ofr: 

H^Hoo < [~(m + 4)] 

Remark. For quadratic splines, one can prove that HS2II00 — 2.5 for all partitions. For 
uniform partitions, one gets Halloo = §[jf ~ 1.4734 

4.2. Near-best dQIs. Let us consider the family of dQIs of degree m defined by 

Qf = Q P J = J2^ B *- 

depending on the two integer parameters p > m and q < min(m, 2p). Their coefficient 
functionals depend on 2p + 1 parameters 

***(/) = E A *(*)m+ S ), 

s=—p 

and we impose that Q is exact on the space U q . The latter condition is equivalent to 
Qe r = e r for all monomials of degrees < r < q. It implies that for all indices i, the 
parameters A;(s) satisfy the system of q + 1 linear equations: 

E X ^ + s = 0<r<q. 

s=—p 

The matrix Vi E M.( g+1 ^ x ( 2p+1 > of this system, with coefficients Vi(r, s) = Q r i+S , is a Vander- 
monde matrix of maximal rank q+ 1, therefore there are 2p — q free parameters. Denoting 
by hi G M. q+1 the vector in the right hand side, with components &j(r) = 0\ , < r < q, 
we consider the sequence of minimization problems, for i £ Z: 

min{||Ai||i; VJAj = fej}. 

We have seen in the introduction that v*(Q) = maxj g z mm ||Aj||i is an upper bound of 
II Qq || 00 which is easier to evaluate than the true norm of the dQI. The objective function 
being convex and the domains being afhne subspaces, we have the following 
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Theorem 2. The above minimization problems have always solutions, which, in general, 
are non unique. 

Let us give an example of optimal dQIs m order to bring out the constraints 

on the partition implied by the optimal property. For m = 2 and p > 2, assume that the 
partition T satisfies, for all i 

Qi-i + 6% < 0%-p + &i+p + di+i, 

then there is a unique optimal solution (here hi = ti — £j_i): 

1 /? 2 1 /? 2 



4 (^i-p — 9i +p )(9i — Oi+p) 4 (8i_ p — 9 i+ p)(9 i+p — 8i) 

K(p) = 1 + 177 .tL 5 v A <*^ = /or s ^ °' 

4 l^i+p — — Vi-p) 

Denote by Q* 2 the associated dQI, then we have [I] 

Theorem 3. The infinite norm of Q* 2 is uniformly bounded independently of p and of 
the partition T: 

IIQp^lloo < 3. 



4.3. Uniformly bounded iQIs. Various types of integral QIs are considered in jT^j [40J [5 . 

Here, we restrict our study to Goodman-Sharma (abbr. GS) type iQIs which appear in 
|21j . They are simpler than those we have already studied in j^U]- Given / = [a, b] and the 
sequence of knots T = {i_ m = ...=to = a <^i<---< t n -i < b — t n — . . . — t n+m }, let 
Bi be the B-spline of degree m and support [tj_ m , t i+ i] normalized by X^o™" 1 = 1, and 
let Mj_i(t) be the B-spline of degree m — 2 with support Ej_i = tj], normalized 

by ffl = Ai( e o) = / Mi-i(t) = 1- The original GS-type iQI can be written as follows 

ra+ m— 1 
i=0 

where //„(/) = /(*o), K+m-i{f) = f(t n ) and, for 1 < i < n + m - 2, 

ik(f) = [ Mi^{t)f{t)dt. 

J a 

It is easy to verify that G\ is exact on ITi and that ||Gi||oo = 1- It has also interesting shape 
preserving properties: for example it is proved in [21] that Gi preserves the positivity and 
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the convexity of the approximated function /. Moreover, the authors also prove that 

G x e 2 -e 2 = (iSiea - e 2 ). 

m + 1 

Let us also prove that it preserves the monotonicity of /. 

Theorem 4. If f is a monotone function on I, then G\f is also monotone with the same 
sense of variation. 

Proof. Let us use the simplified notation G\f — ^S)" 1- PkBi- Then the first derivative 
is DG\f = m y~]7--T~ 1 T—r^ Bf, where B* is the B-spline of degree m — 1, normalized as 
Bi. Assuming that / is monotone increasing, we have to prove that DG\f > 0. For that, 
it is sufficient to prove that the sequence {(//j — fii-x), l<i<n + m — 1} is monotone 
increasing. For i = 1, we have fix — /to = Mof — f(to); as M = 1, the mean value 
theorem gives M f = /(£ ) for some £ e [t ,*i], hence fM 1 -fk = /(£ ) - /(*o) > 0. 
The same kind of proof holds for i — n + m — 1. Now, for 2 < i < n + m — 2, we know 
that Mi - Mi-x = -DB*_ 2 = (see e.g. [7j), therefore & - /Vi = 7-37— J(-DB*_ 2 f) 
and, after integration by parts, ttj — £tj_i = 7-77 — f B*_ 2 f f > 0. 

Therefore, we can see that the operator Gx is very close to the Schoenberg's operator. 
^From the expression of Gxe 2 — e 2 , we can deduce the family of GS-type iQIs defined by 

n+m— 2 

G 2 f = f(t )B + ^2 \ a iP-i-i(f) + hfiiif) + Cifi i+ x{f)}Bi + f(t n )B n+m _x, 

i=X 

which are exact on LT2. The three constraints G 2 ek = e^, k — 0, 1, 2, lead to the following 
system of equations, for l<i<n + m — 2: 

cii + bi + Ci = 1, Qi-xai + 6ibi + 6 i+1 Ci = 0;, fif}x a i + A«i + /4+i c i = 9? '■ 
This is a consequence of the following facts 

/ii(ei) = / tMi-x{t)dt = — y^U_ m+s = 9 h 

2 

ti 2) = ^(e 2 ) = / t 2 M^x(t)dt = — — rr^CZi) 



Kr<s<m 



m(m + 1) 

For m = 2 (quadratic splines), here are the explicit expressions of the coefficients: 



(h^x + hi){hi-.x + K + h i+1 ) ' + h + h i+1 )(hi + h i+1 ) ' 
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and hi = 1 — a,i — Ci. It is clear that HG^Hoo < 1 + 2max(|aj| + |q|) < 5. In fact, this 
upper bound is valid for any degree m: the computation of coefficients and the proof of 
this interesting result will be given in j^j. 

Theorem 5. The iQIs G2 are uniformly bounded independently of the partition and of 
the degree m. More specifically, one has 



5. BlVARIATE QUADRATIC SPLINE QIS ON NON-UNIFORM PARTITIONS 

At the author's knowledge, the only bivariate box-splines which have been extended 
to non uniform partitions of the plane are C 1 -quadratic box-splines on criss-cross trian- 
gulations 13J[3Q.|. Recently jSEj-jSH] we have constructed a set of B-splines generating 
the space of quadratic splines on a rectangular domain and having their support in this 
domain. Moreover, we have defined a discrete quasi-interpolant which is exact on II2 
and uniformly bounded independently of the partition (it is different from the operator 
introduced in [T3|, see e.g. [TOJ , [TT| ) . For the sake of simplicity, we assume here that the 
domain is the whole plane endowed with a nonuniform criss-cross triangulation obtained 
by drawing diagonals in each rectangle R^ = [xj-i, Xj\ x [yj-i, yj] of the partition defined 
by the two sequences X = {xi, i £ Z} and Y = {yj, j G Z}. We set hi = Xi — 

= yj ~~ Vj-ii s i = \( x i-i + x i) an d tj = \{yj-\ + yj)- Let Bij be the B-spline whose 
octagonal support is centered at the point Uij = (si,tj). Let ITj,,' be the continu- 
ous piecewise affine pyramid (of egyptian type) satisfying Yiijiojij) = 1, whose support 
is the central rectangle Rij of £jj. As f R Uij = \h. t kj, we can also define the pyramid 

v = tHkT^v norma hzed by J R Uij = 1. Finally, in the same way, let \ij be the charac- 
teristic function of Rij and let Xij = j^Xij normalized by j R Xij = 1- F° r monomials, 
we use the notation e rs (x, y) = x r y s . We know study some families of discrete and integral 
quasi-interpolants. Detailes proofs will be given elsewhere. 

5.1. Discrete quasi-interpolants. The simplest dQI is the Schoenberg operator defined 



G 



2 00 — 



< 5 



by 
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It is well known ([H], [IE]) that Sie rs = e rs for < r, s < 1 and 

Sie 2 o = e 20 + ^ £ h i B v> Sie ° 2 = e ° 2 + 7 £ ^-Sij- 

(i,i)ez 2 (»,i)ez2 

Si is clearly a positive operator and it preserves the bimonotonicity and the biconvexity 
of / (i.e. the monotonicity and the convexity in the directions of coordinate axes). This 
can be proved by using the expressions of partial derivatives diB^ and d^B^ which are 
piecewise affine functions whose values are given in the technical report . 



5.2. Integral quasi-interpolants. We now study the two following iQIs: 

T l f= 

{i,j)ei? 

where (/, g) = J R2 f(v)g(v)dv, and the GS-type iQI: 

(i,j)6Z 2 

Let ^ij(f) = j^- J R . f be the mean value of / on Rij, then one can also write 



Gif= £ IMjiDBy. 

These two operators are very close to each other and to the Schoenberg operator defined 
before. In particular, they are exact on bilinear polynomials and they satisfy respectively 

7ie 20 = e 20 + ^ Y h i B v> Tie o 2 = e ° 2 + Jo k j B v> 
(i,j)ez 2 {i,i)& 2 

Gie 20 = e 20 + - Y h i B vi G i e 02 = e 02 + ^ £ fcf-By- 

Moreover, as Si, they both preserve the bimonotonicity and the biconvexity of /. 

^From the properties of T\ and G\ on monomials e 2 o and eo 2 , one can deduce the two 
following iQIs which are both exact on the space Il 2 of bivariate quadratic polynomials: 

T 2 f= Y U,M l3 )B l3 , G 2 f= Y l J->-,j)lh r 

(i,j)ez 2 (i,j)& 2 

where the two functions My and ipij are respectively defined by 

Mij = a,iHi-ij + ajiij+ij + bijUij + CjHij-i + cjHij+i 
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The coefficients of these functions are the following 

-3/i? _ -3/i? 



(hi-i + hi)(3hi-i + Ahi + 3/i i+ i) ' (3/ii_i + 4/ij + 3h i+1 )(hi + h i+ i) 
3 kj _ 3 kj 



J (fcj-i + kj) (3^-1 + 4fcj + 3fcj+i) ' J (3fcj-_i + + 3% + i) + ' 
and fey = 1 — (a^ + a, + Cj + Cj), for the operator T%. It is easy to verify that for all 
criss-cross triangulations, one has |aj|, |aj|, \cj\, \cj\ < |. For the operator G 2 , we get 

-/i? -/i? 



+ + + ' (hi-i + hi + h i+ i)(hi + h i+ i) 



-k 2 — k 2 

j — j 



(kj-i + kj){kj-\ + kj + kj+i) ' (kj-i + kj + kj +1 )(kj + k j+1 ) ' 
and fiij = 1 — (otj + + 7j + 7^). It is easy to verify that for all criss-cross triangulations, 
one has \c\i\, |7j|, |tj| < 1. /.From these inequalities and the fact that 

Halloo < 1 + 2max{|aj|, |a f |, |c,-|, |c,-|}, HG2IU < 1 + 2max{|a; i |, |a;|, 

one can deduce the following interesting result. 

Theorem 6. For all non-uniform criss-cross triangulations, the operators T 2 and G2 are 
uniformly bounded. Moreover, 

Halloo < 7 and ||G 2 ||oo < 9. 

Remark. In the case of uniform criss-cross triangulations, one gets respectively: 

— 3 8 _ 11 

a>i — a-i — Cj — Cj — ~2Q-> hj — g> =^ || -Z2 1| 00 ^ "g - ~ 2.22 

—15 7 

a» = ai = 7j = 7j = — , Pij = -, \\G 2 \\oo < - ~ 2.33. 

5.3. Powell-Sabin quasi-interpolants. Recently (|25j.|26j). using an interesting result 
by Dierckx ^Hj, we have introduced and studied new families of quadratic splines quasi- 
interpolants defined on Powell-Sabin type triangulations. As for the previous QIs of this 
section, we have obtained some families of QIs which are both exact on II2 and uniformly 
bounded independently of the partition. 



6. Some applications 
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6.1. Approximation of functions. /From a classical result in approximation theory 
(see e.g. [20J) we know that if Q is an operator defined on a space of smooth functions / 
with values in a space of splines S, one has 

11/ - Q/iloo < (1 + ||Q||oo)rfoo(/, «S) 

As the various QIs studied above have uniformly bounded norms, their approximation 
order is only governed by doo(f,S), i.e. by the distance of / to the maximal space of 
polynomials included in S. Since the values of ||Q||oo are small, we obtain quite good 
approximants which can be used in various fields of numerical anlysis. 

6.2. Approximation of zeros. Quadratic dQIs are simple and good approximants and 
their zeros are rather easy to compute, so they give good approximations of the zeros of the 
approximated function. We already did some computations with orthogonal polynomials, 
and the first results are encouraging (see e.g. [SB] > |SH|)- 

6.3. Quadrature formulas. For the same reasons, quadrature formulas (QF) are easily 
obtained by integrating spline QIs. An interesting univariate example is given in [36J. The 
study of bivariate and trivariate QF is still in progress and the results already obtained 
in two and three variables are also encouraging |T%] . 

6.4. Pseudo-spectral methods associated with quasi-interpolants. Derivatives of 
QIs give quite good approximations of derivatives of the approximated function. This 
simple fact is the basic idea for developing pseudo-spectral methods based on univariat 
and multivariate dQIs with low degrees. 

Keywords: spline operators, quasi-interpolants. 
AMS Classification: 41A36. 
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